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Abstract 

The suucture of a single layer of graphene on Ru(0001) has been studied using surface x-ray diffraction. 
A surprising superstructure has been determined, whereby 25 x 25 graphene unit cells lie on 23 x 23 unit 
cells of Ru. Each supercell contains 2x2 crystallographically inequivalent subcells caused by corrugation. 
Strong intensity oscillations in the superstructure rods demonstrate that the Ru substrate is also significantly 
corrugated down to several monolayers, and that the bonding between graphene and Ru is strong and cannot 
be caused by van der Waals bonds. Charge transfer from the Ru substrate to the graphene expands and 
weakens the C-C bonds, which helps accommodate the in-plane tensile stress. The elucidation of this 
superstructure provides important information in the potential application of graphene as a template for 
nanocluster arrays. 
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The detailed structure determination of single-layer graphene on well-defined surfaces is a sig- 
nificant goal in materials science and solid-state physics - it is most probable that future electronic 
devices based on graphene layers will be fabricated on crystalline substrates 111, hence a knowl- 
edge of how substrates affect graphene is of paramount importance if the latter's structural and 
electronic properties are to be tailored [2]. In addition, it has been recently discovered [13j,|4j,|5j,l6j,LZ0 
that, when grown on crystalline transition metal surfaces, graphene can form superstructures re- 
sulting from moire superpositions of (m x m) carbon hexagons on (n x n) metal surface cells. 
It is still disputed as to whether observed features within these supercells are caused by electron 
density fluctuations over a basically flat structure [6], or whether there is an actual buckling of the 
graphene sheet ||3|,l8|]. A structural clarification would identify the potential of graphene in applica- 
tions such as molecular recognition, single-mo 



for biological or catalytic applications 
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ecule sensing [9], and nanocluster array templates 
[|. 1 1211 . Here we show, using surface x-ray diffrac- 
tion (SXRD), that graphene forms a surprising superstructure when grown on Ru(0001), whereby 
25 x 25 graphene unit cells lie commensurately on 23 x 23 unit cells of Ru. Characteristic 
intensity oscillations in the SXRD data prove that not only the graphene but also the Ru down to 
several atomic layers are significantly corrugated, indicating that the bonding between the single 
graphene layer and Ru is unusually strong. 

The structure of graphene on Ru(0001) has already been investigated using scanning tunneling 
microscopy (STM), conventional electron microscopy, x-ray photoelectron spectroscopy, Raman 
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spectroscopy, and low-energy electron diffraction (LEED) and microscopy |3, 5, fj[ 1, 
well as density functional theory (DFT) Q8Q. Until now it has remained a contentious issue as to 
what registry exists between graphene and the underlying Ru substrate. It has been suggested that 
12 x 12) graphene hexagons sit on (1 1 x 1 1) Ru surface nets (described henceforth as 12-on-l 1) 
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13J], while an ll-on-10 structure has also been proposed [6]. If one assumes an in-plane 



lattice constant for graphene equal to that for graphite (2.4612 A), the former structure would have 
an in-plane tensile strain on Ru (a = 2.706 A) of 0.78 %, while the latter would be compressively 
strained by only 0.05 %. 

The positions of the first order diffraction signals associated with these two superstructures are 
1.100 and 1.0909 in-plane reciprocal lattice units (r.l.u.) of the underlying Ru lattice, i.e., they lie 
within less than 0.01 r.l.u. of one another. With a conventional LEED system, one can only achieve 
accuracies of about one to two percent, excluding an unambiguous identification of the structure 



using this method. Only surface x-ray diffraction is capable of achieving this goal 
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FIG. 1: (color online) (a) A LEED image of the graphene/Ru(0001) surface taken at an electron energy of 
74 eV. (b) An STM image of graphene on Ru(0001), highlighting the supercell containing four subcells. 



Two samples of graphene were prepared on separate occasions on the same sputtered and an- 
nealed Ru(0001) single crystal. In both cases, the crystal was heated to 1115 K in ultra-high 
vacuum (UHV) and a single layer of graphene was deposited by dosing ethene at a pressure of 
2 x 10 5 Pa for three minutes II 1711 . The temperature was then held at 1 1 15 K for a further 60 sec- 
onds. For the first sample, the crystal was then cooled at a rate of 0.4 Ks _1 down to 915 K, and 
from there at a rate of 0.8 K s -1 down to 610 K, after which the heating was turned off. In the case 
of the second sample, cooling was approximately four times quicker. It is noted here that these dif- 
ferent cooling rates were chosen to establish whether they affected the final form of the absorbate 
structure. No significant difference between the two samples could be detected, however. Typical 
LEED and STM images after cooling to room temperature are shown in Fig.Q] 

The samples were transferred in UHV to a minichamber (10~ 7 Pa) equipped with a hemispher- 
ical Be-dome for studies using SXRD [18]. The chamber was mounted on the surface diffractome- 

n 

ter of the Materials Science beamline, Swiss Light Source uM- Structure factors were recorded 
using the Pilatus 100k pixel detector 120(1 . The photon energy was 12.4 keV and the transverse 
and longitudinal coherence lengths were both 1 /urn. 

A calibration of reciprocal space was achieved to two parts in 10000 by using bulk Ru Bragg 
reflections as reference points. A typical set of in-plane scans is shown in Fig. [2] In Fig. (2b), 
satellites on either side of the Ru (01) crystal truncation rod (CTR) can be seen. The primary 
graphene (01) signal [see Fig. |2£d)] at k = 1.087 r.l.u. is, in itself, no proof of a commensurate 
reconstruction, as the graphene could in principle lie incommensurately above the Ru substrate. 
However, the presence also of a signal an equal distance on the other side of the Ru CTR indicates 
that a reconstruction must exist [Fig. |2£c)]. We found several other signals proving a true recon- 
struction elsewhere in reciprocal space [see Fig. Oa)] Note that the (0001) surface of hexagonal 
close-packed systems commonly display sixfold symmetry, although the symmetry of a perfect 
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FIG. 2: (color online) Summary of the diffraction data for the graphene/Ru(0001) system, (a) Schematic 
reciprocal space map showing where data was recorded. The red circular dots indicate points recorded in 
plane at I = 0.4 r.l.u. The in-plane scan along the /c-direction in the neighborhood of the (01) CTR of Ru 
at / = 0.4 r.l.u. shown in (b) is indicated by the green (gray) line. The positions of superstructure rods 
shown in (e) are indicated by the orange diamond and blue square, (c) and (d) High-resolution scans across 
the superstructure signal of graphene on Ru, detailing the 25-on-23 reconstruction. The data were fit to a 
pseudo-Voigt profile (solid black curves), while the positions where 13-on-12 and 12-on-ll reconstruction 
diffraction signals would lie (dotted lines) are also shown. 

hcp(0001) surface is threefold. This apparent increase in symmetry is caused by surfaces con- 
taining regions separated by atomic steps of half a unit cell height, resulting in a 180° rotation of 
adjacent terraces. This is why only one in-plane axis is shown in Fig. [2fa). 

The positions of the two reconstruction signals shown in Fig. (2c) and (d) and those of all the 
other superstructure signals investigated, indicate however, that the superstructure complies with 
neither of those proposed so far. In fact, the signals sit exactly at 21/23 and 25/23 r.l.u., to within 
0.0002 r.l.u. from which it is unambiguously clear from our SXRD data that the reconstruction 
is in fact 25-on-23, that is 25 x 25 graphene honeycombs sitting commensurately on 23 x 23 
Ru unit cells. This signal cannot be explained as originating from the incoherent addition of 
diffraction signals from large domains of 13-on-12 and 12-on-ll supercells, as the linewidth of 
the superstructure signal is, at (5.4 ±0.1) x 10~ 3 r.l.u., significantly narrower than the separation 
of any independent 13-on-12and 12-on-ll signals of 7.6 x 10 3 r.l.u. [see Fig.[2£c) and (d)]. Also, 
a model consisting of a random distribution of 13-on-12 and 12-on-l 1 supercells differs from that 
of the 25-on-23 structure by maximum in-plane displacements of the carbon atoms of less than 
0. 1 A, which are significantly smaller than typical vibrational amplitudes at room temperature, 



4 



hence can be ignored. 

The large extent of this supercell, covering over 33 nm 2 and containing 1250 carbon atoms 
is surprising, since the structure forms in a process involving temperatures above 1000 K, where 
more than 300 eV thermal vibrational energy is stored in the supercell. In the scanning tunneling 
microscope image of graphene on Ru(0001) shown in Fig. QIb), a supercell is highlighted. This 
contains not one, but four parallelogram structures. It is still disputed whether the hill-like features 
are formed by a physical corrugation of the graphene sheet, or is caused by electron density waves 
in an essentially flat graphene layer [6]. Although this cannot be decided from the STM data 
alone, a recent combined DFT/STM study [|80 revealed that the graphene is indeed significantly 
corrugated when deposited on Ru(0001). It is these features that make this system so interesting 
as a potential nanotemplate. 

The four "subcells" within the supercell cannot map translationally onto one another, as from 
the SXRD data it is clear that the number of unit cells of Ru (23) as well of graphene (25) along 
the edges of the supercell are odd, and one is therefore forced to conclude that the graphene 
supercell must consist of four translationally inequivalent subcells. It is noted that, because of 
the presence of 2 x 2 corrugation periods within each supercell, all superstructure peaks with an 
in-plane distance from the Ru signals of p/23 r.l.u., where p is an odd integer, are systematically 
absent. 

The 21/23 and 25/23 superstructure rods (SSRs) are shown in Fig[2e). Because they only 
provide information on surface reconstructions of the graphene and uppermost Ru-layers, and not 
on bulk properties, we are able to infer important properties of the surface region immediately 
from their qualitative features. First, the signal intensity is strongly modulated, with a periodicity 
of approximately 1.0 r.l.u. (with respect to Ru) in the out-of-plane direction. This modulation can 
only occur if the ruthenium is physically corrugated, that is, if the graphene imposes vertical strain 
[8|]. Importantly, all the maxima have widths of approximately 0.25 r.l.u., which means that the 
Ru substrate must also be significantly corrugated down to about 4 unit cells, or over 1.5 nm. 

Unsurprisingly, the 25/23 rod has significant intensity at low /-values, as at / = 0, this corre- 
sponds to the (010) in-plane graphene peak, which is known to have non-zero intensity. However, 
extrapolation of the 21 /23-rod to I = strongly indicates that it also has non-zero intensity here, 
which can only occur if there are in-plane movements of the atoms within the supercell. 

In order to estimate the corrugation amplitudes and depths of the ruthenium, we used a simple 
model to simulate the two superstructure rods of Fig|2te)- The vertical displacement field of the 
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FIG. 3: (color online) Simple parametric model of the graphene-Ru(OOOl) supercell. (a) The vertical dis- 
placement field of the graphene corrugation. The scale is given in A above the Ru substrate, (b) Simulated 
21/23 and 25/23 superstructure rods of corrugated graphene on Ru(0001), incorporating the qualitative 
features extracted from the experimental data. The graphene has a peak-to-peak corrugation amplitude of 
1.5 A, while that of the uppermost Ru atomic layer is 0.20 A. The corrugation amplitude of the Ru drops 
exponentially by a factor of 0.75 for each successive atomic layer (0.56 per unit cell depth). The minimum 
graphene-Ru distance is 2.2 A, and the mean vertical positions between all the Ru atomic planes assume 
bulk values. 

graphene corrugation [Fig. Oa)] is generated by first calculating a subfield for each of the two 
inequivalent carbon atoms in the conventional graphite unit cell, whereby the vertical distance to 
the Ru-substrate is proportional to the in-plane separation from the nearest Ru-atom. The com- 
ponent fields are interpolated and then added to produce the final displacement map. The model 
incorporated only four parameters, namely the graphene corrugation amplitude and the minimum 
graphene-Ru distance, for which we used fixed values determined by DFT calculations [8], and 
the corrugation of the uppermost Ru atomic layer and the exponential decay depth of the Ru- 
corrugation, which were varied to best match the experimental linewidths and intensities. Any 
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in-plane movements of either the graphene or the Ru were ignored, as this would add significant 
complexity to the model, and our a priori knowledge of such movements is very limited. The 
simulation is shown in Fig. [3] Despite the simplicity of this model, the qualitative agreement is 
impressive. The peak-to-peak corrugation amplitude of the uppermost ruthenium atomic layer 
is 0.20 A, which decays exponentially with depth, with a characteristic length of 1.7 unit cells 
(3.4 atomic layers). Note that the graphene and ruthenium corrugations were chosen to be in 
phase (i.e., peak above peak and valley above valley). If the corrugations are made to be out of 
phase (peak above valley, valley above peak), the agreement with the experimental data is poorer. 

The effect on the shape of the SSRs of the graphene corrugation amplitude is fairly insensi- 
tive, due to the low x-ray scattering amplitude of carbon. I(V) curves from low-energy electron 
diffraction may provide important further quantitative information regarding the detailed graphene 
structure, due to its higher surface sensitivity. A rigorous fit would also include in-plane move- 
ments of both the carbon and Ru atoms and the possibility that the average height of each Ru 
atomic layer can vary. Allowing in-plane movements could significantly change the magnitudes 
of the calculated corrugations, hence their values given here are still tentative. 

The large depth to which the ruthenium substrate is perturbed is, however, a robust parameter, 
and is indicative of an attendant strong chemical bonding of graphene to the metallic substrate via 
the carbon p z orbitals, which cannot arise through van der Waals bonds. This has recently been 
predicted by DFT calculations, that find strong charge redistributions and a minimum graphene- 
Ru distance of only 2.2 A, incompatible with van-der- Waals interactions [8]. Strong interactions 



have already been seen for graphene on Ni(lll) [|17l. l2l|l . This increased bonding to the substrate 
is connected with expanded and weakened C-C bonds, as indicated by the softened phonons of 
the graphene layer fl 1 3M . 

This in-plane expansion of the graphene layer might therefore accommodate the apparent ten- 
sile strain of the 25-on-23 supercell when assuming a bulk-like graphite in-plane lattice constant 
for graphene. Indeed, allowing for an expansion of the C-C bond of approximately 1 % due to 
electron transfer would result in nominally zero heteroepitaxial strain for this 25-on-23 structure, 
although, of course, because we observe elastic deformation of the Ru substrate, there must be 
stress in the graphene layer that imposes strain both in itself and in the substrate. In other words, 
the Ru-C bonding causes the graphene to dilate in plane and the 25-on-23 structure to have the 
lowest surface energy. 

In conclusion, the structure of the supercell of graphene on Ru(0001) has been elucidated and 
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been shown to consist of 25 x 25 unit cells of graphene on 23 x 23 unit cells of Ru. On the one 
hand, this large in-plane extent of over 6 nm, and on the other, the large corrugation amplitudes 
of the Ru-substrate, which indicate a strong bond between the graphene and ruthenium, suggests 
that this system may be ideally suited as a robust template for arrays of nanoclusters or macro- 
molecules. 
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